Why could Electron Spin Resonance be observed in a heavy fermion Kondo lattice? 



B.I. Kochelaev, 1 S.I. Belov, 1 A.M. Skvortsova, 1 A.S. Kutuzov, 1 
J. Sichelschmidt, 2 J. Wykhoff, 2 C. Geibel, 2 and F. Steglich 2 

1 Theoretical Physics Departement, Kazan State University, 420008 Kazan, Russia 
2 Max Planck Institute for Chemical Physics of Solids, 01187 Dresden, Germany 

(Dated: July 22, 2009) 

We develop a theoretical basis for understanding the spin relaxation processes in Kondo lattice 
systems with heavy fermions as experimentally observed by electron spin resonance (ESR). The 
Kondo effect leads to a common energy scale that regulates a logarithmic divergence of different spin 
kinetic coefficients and supports a collective spin motion of the Kondo ions with conduction electrons. 
We find that the relaxation rate of a collective spin mode is greatly reduced due to a mutual 
cancelation of all the divergent contributions even in the case of the strongly anisotropic Kondo 
interaction. The contribution to the ESR linewidth caused by the local magnetic field distribution 
is subject to motional narrowing supported by ferromagnetic correlations. The developed theoretical 
model successfully explains the ESR data of YbRli2Si2 in terms of their dependence on temperature 
and magnetic field. 

PACS numbers: 76.30.Kg, 71.27.+a, 75.40Gb 



The discovery of electron spin resonance (ESR) in the 
heavy fermion Kondo lattice YbRh 2 Si2 at temperatures 
below the thermodynamically determined Kondo tem- 
perature (Xk = 25 K) became a great surprise for the 
condensed matter physics community H According to the 
common belief, based on the single ion Kondo effect, the 
ESR signal should not be observable by at least two rea- 
sons. Firstly, at temperatures below Tk the magnetic 
moments of the Kondo ions should be screened by the 
conduction electrons; secondly, the ESR linewidth was 
expected to be of the order of Tk = 500 GHz. The ex- 
perimental results were completely opposite: at X-band 
(9.4 GHz) and T = 0.7 K a linewidth of 0.3 GHz was 
observed and the ESR intensity corresponds to the par- 
ticipation of all Kondo ions with a temperature depen- 
dence following a Curie- Weiss law!^ Moreover, the angu- 
lar dependence of the resonant magnetic field reflects the 
tetragonal symmetry of the electric crystal field at the 
position of the Yb-ion with an extremely anisotropic g- 
factor (g± = 3.56, g\\ = 0.17 at T = 5 K). Similar results 
were obtained later for Ybh^Si^ and the main features 
of the ESR phenomenon were confirmed recently at very 
high frequencies up to 360 GHzSl Up to now this para- 
dox was not resolved on the basis of the properties of 
these materials studied by other methods. In particu- 
lar, the heavy- fermion system YbRli2Si2 attracted a lot 
of attention due to the existence of an antiferromagnetic 
quantum critical point (QCP), at which antiferromag- 
netic order disappears smoothly as T — > 0, by varia- 
tions of an external magnetic field, pressure or doping. 
Near a QCP and up to surprisingly high temperatures, 
a new phase appears exhibiting non-Fermi liquid (NFL) 
behavior. For YbRli2Si2 the electrical resistivity linearly 
increases with temperature and the Sommerfeld coeffi- 
cient of the electronic specific heat diverges logarithmi- 
cally upon cooling down to 0.3 K. At lower temperatures 
this divergence becomes even stronger. It seems that 
many properties of heavy fermions in the NFL state of 



YbRli2Si2 can be described in terms of quasi-localized /- 
electrons. In particular, the underlying fluctuations near 
the QCP can be considered as locally critical, having an 
atomic length scale^. On the basis of a localized 4/ 
electron approach we successfully studied the tempera- 
ture dependence of the static magnetic susceptibility of 
YbRli2Si2 and Yblx2Si2 at temperatures below Tj£Q 

In this paper we show that the key ingredient of the 
ESR signal existence in a heavy fermion Kondo lattice in 
the NFL state is a formation of a collective spin mode 
of quasi-localized /-electrons and wide- band conduction 
electrons even in a strongly anisotropic system. We shall 
discuss, also, the role of the local ferromagnetic fluctua- 
tions in YbRhaSiiP and YbIr 2 Si 2 il . 

Our basic theoretical model includes the kinetic en- 
ergy of conduction electrons, the Zeeman energy, the 
Kondo interaction of the Yb 3+ ions with conduction elec- 
trons, and the coupling between the Yb 3+ ions via con- 
duction electrons (the RKKY interaction). The lowest 
multiplet of the free Yb 3+ ion is 2 F 7 / 2 with total mo- 
mentum J — 7/2. The tetragonal crystal electric field 
splits this multiplet into four Kramers doublets that are 
well separated from the lowest (ground) doubletP^ such 
that the physics of low energy spin excitations can be de- 
scribed by the lowest doublet. After projection onto the 
Kramers ground state we obtain an effective Hamiltonian 
TL = Hq + TLsa + Hrkky with 
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Here a labels the orientation of the conduction electron 
spin, Sj is the spin- 1/2 operator of the j-th Yb-ion, B is 
the external magnetic field multiplied by the Bohr mag- 
neton; <7|i , g± are the (/-factors for B aligned parallel and 
perpendicular to the crystal symmetry axis; J and Iij are 
the Kondo- and RKKY-coupling constants; g a and (Tj 
denote the g-factor and the spin operator of conduction 
electrons at the position of the Yb-ion. In the following 
we shall discuss the contributions of the magnetic dipole- 
dipole and the spin-phonon interactions and the role of 
the translational diffusion of the /-electrons. 

As a first step we find the renormalized relaxation rate 
of the transverse (relative to the external static magnetic 
field) magnetic moment of the Kondo ions toward the 
conduction electrons which remain in a thermodynami- 
cal equilibrium state (this is the Korringa relaxation rate 
T ss ). For an isotropic system with g± = gu = g the re- 
sult in second order in TL sa is well known: in terms of 
Eq. T ss = I (gpJ) 2 k-sT . However, in the case of 
an antiferromagnetic coupling, J > 0, second order at 
low temperatures is not sufficient, and we have to im- 
prove perturbative calculations. This can be done by the 
renormalization of the coupling constant in the spirit of 
the Andersons "Poor Man's" scaling.^ We provide our 
results for the static magnetic field oriented perpendic- 
ular to the crystal symmetry axis only, since the most 
detailed experiments with YbRli2Si2 were performed for 
this orientation and, moreover, gu data of YbRli2Si2 were 
not accessible with the available experimental equipment 
due to gu < 0.2. The renormalized Korringa relaxation 
rate was found for k B T > B using the method of func- 
tional variational derivatives: 
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The parameter tp reflects the logarithmic behaviour at 
temperatures close to the characteristic value Tqk- For 
the Kramers ground state we have found 
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with W as the band width of the conduction electrons. 
One can see that in the isotropic case the standard result 



T K = Wexp 



can be obtained asymtotically for 
(j5|. As expected, the Korringa re- 
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laxation rate is logarithmically divergent upon lowering 
the temperature to T GK ((p -> 0) : T ss cx 1/ln 2 (T/T GK ). 
The Overhauser relaxation rate T aa (magnetic moment 
of the conduction electrons relaxes toward the Yb 3+ spin 
system, being in the equilibrium with a thermal bath) 
can be found from the relation 

T ss /T aa = 2pk B (T + e±)g a /g± . 0J- = E 7 «- ( 6 ) 
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Here 9± is the Weiss temperature, which originates from 
the RKKY interaction in the molecular field approx- 
imation. Evidently, T aa is also divergent as T aa oc 
1/ In 2 (T/Tqk)- At first glance, these results confirm 
the commonly accepted belief that the ESR linewidth of 
Kondo ions (as well as conduction electrons) is expected 
to be too large for its detection. However, an equilibrium- 
state approximation for the conduction-electron-spin sys- 
tem is not valid to study the ESR response of the samples 
with a high concentration of Kondo ions^ which is espe- 
cially the case for the Kondo lattice. Instead one has to 
derive coupled kinetic equations for both magnetic mo- 
ments of Kondo ions and conduction electrons. 

The kinetic equations of motion for the transverse mag- 
netizations of localized moments and conduction elec- 
trons are coupled by two additional kinetic coefficients 
r s(J and r CTS , respectively. For a correct analysis of a 
stationary solution one has to introduce, also, relaxation 
rates of the Kondo spin system and conduction electrons 
toward the thermal bath ("lattice") T s l and T^l, re- 
spectively. For the anisotropic system such equations 
were derived in second order of the Kondo interaction by 
the methods of nonequilibrium statistical operator and 
Green functions^. The renormalized kinetic coefficients 
r scr and T as were found for k B T > B to be: 
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r s(T can be found from a relation similar to Eq. ([6]) . One 
sees that both coefficients are also divergent in the same 
way (via the parameter tp) as the Korringa and Over- 
hauser relaxation rates. To study the ESR response of 
the total system, one has to look for the poles of a solu- 
tion of coupled equations after their time Fourier trans- 
form. As a result we obtain two complex frequencies: 
their real parts give resonant frequencies, their imaginary 
parts give the corresponding relaxation rates. We are in- 
terested in a solution close to the Kondo-ion resonance 
frequency. 

The coupling between two systems is especially impor- 
tant if the relaxation rate of conduction electrons toward 
the Kondo ions is much faster than to the lattice and 
the resonant frequencies are close to one another ("bot- 
tleneck" regime): T arT ^> T^l, \^>s — w<t\- It is well known 
that in the case of an isotropic system and equal Larmor 
frequencies, the ESR linewidth in the bottleneck regime 
is greatly narrowed due to conservation of the total mag- 
netic moment (its operator commutes with the isotropic 
Kondo interaction and the latter disappears from the ef- 
fective relaxation rate). The same situation remains, if 
one takes into account the Kondo effect^. In the op- 
posite case of a strongly anisotropic Kondo interaction 
the results obtained in second order do not show any 
sufficient narrowing of the ESR linewidth in the bottle- 
neck regime^. However, the renormalization of all ki- 
netic coefficients makes the situation completely differ- 
ent. The relaxation rate of the collective mode with a 
frequency close to the Kondo-ion resonance now follows 
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for T > T GK : 

r C on = r sL + rf + rf L + b 2 f(t) ; (8) 
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The simplified expressions for an effective Korringa re- 
laxation rate and an effective relaxation rate of con- 
duction electrons to the lattice in the case <p <1 are: 
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B 2 F(T) describes the external magnetic field dependence 
which is a consequence of a partial opening of the bot- 
tleneck at large magnetic fields due to different g-factors 
of the Kondo ions and conduction electrons. The co- 
efficient F(T) is a rather extended expression contain- 
ing the temperature dependence, ft is important to note 
that, instead of being divergent, the effective Korringa 
relaxation rate r^f is greatly reduced and goes to zero 
at T — > Tqk ! This is a remarkable result: although 
the total magnetic moment does not commute with the 
strongly anisotropic Kondo interaction at all, the diver- 
gent parts of all the kinetic coefficients cancel each other 
in the collective spin mode due to the existence of the 
common energy scale Tqki regulating their temperature 
dependence at T — > Tqk ■ The effective relaxation rate of 
conduction electrons to the lattice ^ is also greatly re- 
duced, becoming proportional to temperature and mim- 
icking the usual Korringa relaxation rate. These results 
allow one to conclude that the main reason of the ob- 
servability of ESR in a Kondo lattice is a formation of a 
collective spin mode in the bottleneck regime and in the 
presence of a Kondo effect. Another important ingredi- 
ent - the short-range ferromagnetic fluctuations due to 
the RKKY interactions will be discussed below. 

Now we have to consider the broadening of the ESR 
linewidth which is represented by the kinetic coefficient 
r s L- An obvious contribution comes from the distribu- 
tion of effective local magnetic fields due to spin-spin in- 
teractions of the Yb-ions and a variation of the ^-factors 
due to defects of the lattice. In particular, the usual 
magnetic dipole-dipole interactions yield approximately 
Ai/joc ~ 700 Oe, while the observed ESR linewidth in 
YbRh 2 Si 2 at the X and Q bands is AH ESR ~ 200 Oe at 
T = 5 K. The contribution from the RKKY interactions, 
which become highly anisotropic after projection onto 
the Kramers ground state, is expected to be much larger. 
Therefore, it is evident that some narrowing mechanism 
for these type of contributions should exist. 

ft is well known that the broadening of the ESR line 
by the distribution of local fields can be reduced in the 
bottleneck regime due to fast reorientation of the Kondo 
ion moment caused by the Korringa relaxation^. How- 
ever, we expect that in the Kondo lattice the motional 
narrowing due to translational diffusion of quasilocalized 
/-electrons in the NFL state can be much more effec- 
tive. An elementary step of the diffusion is a jump of 



an /-electron (or a hole in the closed /-shell) from one 
Yb site to the next nearest site. At the same time this 
jump can be easily done, if the local field, created by the 
RKKY interaction, is the same. This can happen only, 
if the RKKY interaction for the nearest Yb ions is fer- 
romagnetic, ft is worth to mention that the importance 
of ferromagnetic fluctuations for the ESR observability 
at T < Tk in Kondo lattice with he avy fe rmions was 
discussed recently in a different contextP^. 

An additional broadening of the ESR linewidth can 
appear due to the spin-phonon interaction of the Kondo 
ions. The main contribution comes from the two-phonon 
Raman and Orbach processes at temperatures above a 
few K. The Orbach process via the excited energy level 
at a given crystal field splitting A yields the temperature 
dependence r° L rbach cx [exp(A/fc B T) - f]" 1 . The same 
temperature dependence is provided by the scattering of 
conduction electrons via the excited energy level. 

Concerning the ESR resonance frequency of the col- 
lective mode, the situation is somewhat different. For a 
single Kondo ion it is well known that besides the usual 
Knight shift of the ESR resonance frequency, the Kondo 
effect results in a divergent logarithmic term. The same 
happens with the resonance frequency of the conduction 
electrons. We have found that all divergent parts of the 
ESR resonance frequency cancel each other in the collec- 
tive mode similar to the relaxation rate described above. 
However, the RKKY interaction provides an additional 
local field at the Yb-ion and the Weiss constant 6± in the 
spin susceptibility. In the molecular field approximation 
involving both the Kondo and RKKY interactions, 6± 
becomes also subject of the Kondo renormalization. As 
a result, the ESR resonance frequency contains the diver- 
gent logarithmic term even for the collective spin mode. 
For the corresponding effective (/-factor g e ^ we obtain the 
following relation (with the magnetic field perpendicular 
to the symmetry axis; less important non-divergent terms 
are omitted): 
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Next, we compare our theory with experimental re- 
sults. At first, we use the divergent logarithmic term in 
g°^ to reveal the characteristic temperature Tqk- Start- 
ing values of (/-factors were taken from the crystal field 
consideration (g± = 3.66)P, the density of states can be 
related to the band width of the conduction electrons 
as p — l/W. The result of the fitting is given in Fig. 
[jj, with pj = 0.05, 6> ± = 0.18 K and T GK = 0.36 K; 
the latter is by two orders of magnitude smaller than 
the Kondo temperature Tk derived thermodynamicalljffS 
and by transport measurements^. Although the revealed 
value Tqk — 0.36 K could be rather approximate, we 
used it to fit the temperature and frequency dependen- 
cies of the ESR linewidth with help of the Eqs. ( 8]9 | in 
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FIG. 1: Temperature dependence of Q-band (a) g-factor and 
(b) ESR relaxation rate T of YbRboS^ (open circles). Solid 
lines denote data fitting: g'f* (Eq. ( |10| ) and rtheor (with two 
contributions as indicated). Inset: Frequency dependence of 
rthcor (solid line) fitted to T(T = 10 K) (data from RefsHQ. 



order to see whether our theory is selfconsistent : 

rthcor = r coU + r s ° L rbach + Const. (11) 

Here Const, represents the local field distribution which 
is greatly reduced by the motional narrowing mechanism 
as discussed above. The results are given in Fig. [T]d. 
The fitting of the temperature dependence of the ESR 
linewidth gave pj = 0.048 and A = 198 K. The latter 
coincides with the first excited energy level of the Yb- 
iorfHJl, confirming that the Orbach processes dominate in 
the spin-phonon relaxation (our estimation of the spin 
relaxation with the conduction electrons via excited level 
gave a significantly smaller contribution), ft can be as- 
sumed that the effective relaxation rate T° ^ of conduc- 
tion electrons to the lattice is negligible. 



From the inset of Fig. [T]d one can see that r t heor in- 
creases roughly as v 2 , in full accord with the experimental 
observation. From the frequency dependence of the ESR 
linewidth we obtain pj = 0.07; different fittings resulted 
in only slight variations of pj . Crude approximations like 
p = 1/W may already explain such variations. Having 
the value pj, we can now estimate the Korringa relax- 
ation rate without the bottleneck regime. According to 
Eq. Q the value pj = 0.05 yields T ss = 51 GHz at 
5 K. This would leave no chance to observe the ESR sig- 
nal, neither at X- or Q-bands nor at higher frequencies 
without formation of the collective spin mode. It is re- 
markable that the Kondo effect, being responsible for a 
suppression of the ESR signal on paramagnetic impuri- 
ties in a metal at low temperatures T < Tk , crucially 
supports it in the Kondo lattice due to the formation of 
the collective spin mode with a dramatic narrowing of the 
ESR linewidth, even in the case of a strongly anisotropic 
Kondo interaction. Recently, additional experimental ar- 
guments in favor of the collective spin mode were given 
in Ref. [201 We should mention that the relation be- 
tween the thermodynamical Kondo temperature Xk and 
our characteristic temperature for the ground Kramers 
doublet Tqk remains an open question in our consider- 
ation. Probably, it could be related to the Kondo reso- 
nance narrowing considered quite recently.^ 

In conclusion, we have revealed new features in the 
properties of the Kondo lattice with heavy fermions and 
have found an explanation of the ESR observability in 
YbRh 2 Si2, that was in the last years a subject of sharp 
discussions in the literature. 
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